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, Aims. We study the origin of large abundances of complex organic molecules in the Galactic center (GC). 
^ . Methods. We carried out a systematic study of the complex organic molecules CH 3 OH, C 2 H 5 OH, (CH 3 ) 2 0, HCOOCH 3 , HCOOH, CH3COOH, 
Oh H 2 CO, and CS toward 40 GC molecular clouds. Using the LTE approximation, we derived the physical properties of GC molecular clouds and 
jQ the abundances of the complex molecules. 

5^ ■ Results. The CH3OH abundance between clouds varies by nearly two orders of magnitude from 2.4xl0~ 8 to l.lxl0~ 6 . The abundance of the 
I other complex organic molecules relative to that of CH3OH is basically independent of the CH3OH abundance, with variations of only a factor 
(73 4-8. The abundances of complex organic molecules in the GC are compared with those measured in hot cores and hot corinos, in which these 
ILJ complex molecules are also abundant. We find that both the abundance and the abundance ratios of the complex molecules relative to CH3OH 
. I in hot cores are similar to those found in the GC clouds. However, hot corinos show different abundance ratios than observed in hot cores and 

■ in GC clouds. The rather constant abundance of all the complex molecules relative to CH3OH suggests that all complex molecules are ejected 
S-h ' from grain mantles by shocks. Frequent (~10 5 years) shocks with velocities >6km s are required to explain the high abundances in gas phase 
5^ °f complex organic molecules in the GC molecular clouds. The rather uniform abundance ratios in the GC clouds and in Galactic hot cores 

indicate a similar average composition of grain mantles in both kinds of regions. The Sickle and the Thermal Radio Arches, affected by UV 
radiation, show different relative abundances in the complex organic molecules due to the differentially photodissociation of these molecules. 
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1. Introduction 



of these molecules from grain mantles. 



Gas-grain interaction in the interstellar medium (ISM) can 
have a large impact on its chemistry because of the desorption 
and/or depletion of molecules onto grains mantles. The gas 
phase abundance of molecules, which are believed to be 
efficiently formed on dust grains such as, methanol (CH3OH), 
ethanol (C2H5OH), and formaldehyde (H2CO), can be en- 
hanced by orders of magnitude due to the ejection/evaporation 
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* Table 7 is only available in electronic at the CDS via 
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One of the interesting features of the gas-phase chem- 
istry after ejection of these molecules is the fast conversion 
of these parent molecules into daughter molecules. The rel- 
ative abundance of dimethyl ether ((CH3)20) is expected to 



time after CH3OH i s ejected to gas phase 



( Millar etal 


Jll99ll 


Horn et alJ 


2004) 



Basically all models of alcohol-driven chemistry have indicated 
that (CH 3 ) 2 and methyl formate (HCOOCH3) will reach their 
largest abundances 10 4 -10 5 years afte r the e jection of CH3OH 
into gas phase. However, iHorn et all d2004l) have shown that 
HCOOCH3 formation in gas phase from CH3OH is less effi- 
cient than previously considered in the models. This casts some 
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doubts on the gas-phase production of HCOOCH3, and they 
propose that this molecule is also formed on grain mantles. 
Th e gas phase/grain fo rmation of (CI^^O has been studied 
by IPeeters et al.l (120061) . showing that the main path to form 
(CH3)20 is by gas-phase reactions and that its formation on 
grains is a minor source of the observed abundances. Other or- 
ganic molecules like formic acid (HCOOH) can be pro duced in 
both gas phase and grain mantle chemistry JLiu et all 2002). 

So far, large abundances of complex molecules have been 
detected in three different ki nds of objects: hot cores associated 
with massive star formation (Ikeda et al. 2001), hot corinos as - 
sociated with low-mass star formation (Bottinelli et al. 2006), 
and t he Galactic center (GC) clouds ( Martin-Pin tado et al.l 
l200ll) . Their physical properties are very different. While the 
hot cores and hot corinos are small (<0.1pc), hot (>100K), 
and very dense (>10 5 cm~ 3 ) condensations, the GC molec- 
ular clouds show averaged scales of 20-30 pc, kinetic tem- 
peratures ~50-200 K, a nd an averaged density of ~10 4 cm" 3 
(Giisten & Philipp 2004). The molecular gas in the inner re- 
gion of our Galax y presents a different chem istry than in 
the Galactic disk fflylartm-Pintado et alJ fl997l) . In the GC, 
widespread large abundances of grain-processed molecules like 
CH3OH, C?Hs OH, and silicon monoxide (SiO) are observed 
jGottlieb et atl 1 19791 iMinh et alJ 1 19921 IMartm-Pin tado et alJ 
Il997l 1200 it lHuttemeisteretal.lll998ir it has been proposed 
that the sputtering of grains and grain mantles produced by 
widespread shocks with moderated velocities of <40 km s 
are the responsible of the "rich" chemistry observed in the 
GC clouds. The origin of the large scale shocks is so far un- 
clear. These shocks could be produced by cloud-cloud colli- 
sions associated wit h the large-scale dynamics in the con text 
of a barred potential (lHasegawa et all 19941 iHiittemeister et alJ 
Il998t iRodnguez -Fernandez et al. 2006), by wind-blown bub- 
bles driven by evolved massive stars JMartin-Pintado et alJ 
I1999I) . or by hydrodynamic (H P) or magn eto hydrodynamic 
turbulence (MHD) jMorris & Serabvnl 19961) . Since the forma- 
tion of some complex organic molecules is believed to pro- 
ceed in gas phase after the passage of the shocks and their 
abundances rapidly evolve, one could use the abundances of 
these molecules to gain insight into the chemistry of complex 
mol ecules in the ISM and to con strain the age of shocks in the 
GC JMartin-Pintado et alJ20ol . 

Following these ideas we made a systematic study of 40 
molecular clouds in the GC region, between Sgr B2 and Sgr C, 
in complex molecules believed to be formed on grains like 
CH3OH, C2H5OH, and HCOOCH3, in gas phase from CH3OH 
like (CH3)20 and in gas phase and/or grains like HCOOH 
and H2CO. We also searched for the less abundant isomer of 
HCOOCH3, the acetic acid (CH3COOH). As tracers of the total 
column density of the molecular gas, we observed molecules 
like C ls O, 13 CO, and CS. Our systematic study shows that all 
the complex organic molecules present similar relative abun- 
dances with respect to CH3OH, except for the regions where 
photodissociation could be important. Furthermore, the abun- 
dance and abundance ratios of these complex molecules in the 
GC are similar to those observed in hot cores, while hot corino 
abundances seem to have similar ratios for different objects. 
Our results suggest that all complex organic molecules have 



Table 1. Molecular line parameters 



Molecule 


Transition 


Frequency 
(MHz) 


E u /k 
(K) 




CH^OH 


3o - > 2o E 


145093.75 


27.06 


2.38 




3_i -» 2_, E 


145097.47 


19.52 


2.11 




3 -> 2 A+ 


145103.23 


13.94 


2.38 




5o -> 4 E 


241700.22 


47.95 


3.91 




5_i -^4_, E 


241767.22 


40.41 


3.75 




5 4 A+ 


241791.43 


34.83 


4.94 


13 CH 3 OH ... 


3(b — * 2o2 E 


141595.48 


26.71 


2.43 




3-i3 — * 2_i2 E 


141597.06 


19.21 


2.16 




3 03 -» 2 02 A+ 


141602.53 


13.60 


2.43 


C 2 H,OH 


4l4 -» 3 03 


90117.61 


9.36 


5.35 




7o7 -> 616 


104487.26 


23.26 


8.61 




9o9 — » 8]g 


142285.05 


37.17 


12.66 


HCOOCH3 .. 


725 -» 6 2 4 E 


90145.69 


19.69 


17.00 




725 -» 6 24 A 


90156.48 


19.67 


17.00 




808 - * 7o7 E 


90227.61 


20.09 


20.99 




808 - * 7o7 A 


90229.63 


20.07 


20.99 




9 45 -> 844 E 


1 1 1408.48 


37.27 


19.13 




9i 8 ^8„E 


111674.10 


28.15 


23.12 




9is -> 817 A 


111682.19 


28.13 


23.12 


(CH3) 2 O fl .... 


726 — » 7 17 


104703.30 


31.07 


5.43 




7o7 -> 6l6 


111783.01 


25.26 


6.80 


HCOOH 


422 — * 321 


90164.25 


23.53 


5.78 




5o5 - » 4o4 


111746.79 


16.13 


9.65 


CH3COOH .. 


8,8 -> 7, 7 A 


90246.25 


20.30 


43.20 




8,8 -> 7, 7 E 


90203.44 


20.30 


43.20 




lO.io -» 9,9 A 


111548.53 


30.50 


54.80 




IO.10 -» 9.9 E 


111507.27 


30.50 


54.80 


CS 


3 ->2 


146969.03 


14.12 


11.56 




5 ->4 


244936.64 


35.28 


19.84 


C 18 


1 -> 


109782.17 


5.28 


0.01 




2 -» 1 


219560.35 


15.82 


0.02 


13 CO 


1 -> 


110201.35 


5.30 


0.01 


H' 3 CO 


2o2 —> loi 


141983.74 


7.11 


10.87 



8 EE, AA, EA, and AE substates blended. Only the most intense 
transition is given. 



been ejected from grain mantles and that their abundances rep- 
resent, in first approximation, the grain mantle composition. 

2. Observations and results 

The observations of all transitions of the complex organic 
molecules shown in Table ^were carried out with the IRAM 
30-m radio telescope at Pico Veleta (Spain). The data were ob- 
tained in different seasons between 1996 and 2003. The half- 
power beam width of the telescope was 24", 17", and 12" for 
the 3, 2, and 1.3 mm bands. The receivers, equipped with SIS 
mixers, were tuned to single sideband with image rejections 
^ 10 dB. The typical system temperatures were 300, 500, and 
900 K for the 3, 2, and 1 .3 mm lines, respectively. We used two 
filterbanks of 256 x 1 MHz and one of 512 x 1 MHz as spec- 
trometers. The velocity resolutions provided by the filter banks 
were 3, 2, and 1.3 km s _1 for the 3, 2, and 1.3 mm bands, re- 
spectively. Spectra were calibrated using the standard dual load 
system. We used the antenna temperature scale (7^) for the line 
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Fig. 1 . Position of the observed sources superimposed on the SiO J=l— >0 maps from the GC region by ( Mar tin-Pintado et alJ 
119971) . 



intensities because the emission is rather extended and com- 
pletely fills the beam. 

Figure ^ shows the location of the sources we ob- 
ser yed superimposed on the la rge-scale SiO map of the GC 
by iMartin-Pintado et al.l 0997I). The sources were selected 
from the SiO maps of IMartin-Pintado et al .] il997h . some 
of these sources had already been observed in C2H5OH by 
Martin-Pintad o et all J200ll) . We also inclu ded some clouds in 
Sgr C , the Thermal Ra dio Arches (TRA^ JSerabvn & Giistenl 
Il987h . and the Sickle dSerabvn & Giistenl Il99lh. where the 
C?HsO H emission was not detected bv IMartin-Pintado et alJ 
(2001). The selected molecular clouds are spread over the re- 
gion between Sgr C and Sgr B2. The positions of the sources 
are given in galactic and equatorial coordinates in Table |2] 
The nomenclature used to d esignate them is similar to that of 
IMartin-Pintado et alJ J200ll) : MC stands for Molecular Cloud 
followed by a G and the galactic coordinates, as recommended 
by the IAU. We grouped the sources in different regions as 
shown in Fig.^and Table [2 

Figures |2] and [5] show samples of line profiles for all 
observed molecules toward selected sources. Only for some 
sources like MC G-0.02-0.07 we did observe all the molec- 
ular lines. The line profiles of molecules such as CH3OH are 
different from those of the CO isotopomers and CS due to the 
blending of several transitions. The complex organic molecules 
were only detected in some of the velocity components ob- 
served in C 18 0, 13 CO, and CS. A good example of this dif- 
ference in line profiles is shown in Figs. |3] and for the source 
MC G+0.76-0.05 where C I8 0, 13 CO, and CS show three or 
more velocity components, but C2H5OH only shows one. Each 
velocity component has been treated as an independent molec- 
ular cloud. Multiple velocity components are also found toward 
MC G+0.13+0.02, MC G+0.18-0.04, and MC G-0.32-0.19. 



3. Analysis of the data 

To derive the physical and chemical properties of the selected 
molecular clouds, we fitted Gaussian profiles to all detected 
molecular lines for each velocity component. The observed 
spectra of some organic molecules are more complex than a 
single Gaussian profile because of the overlap of several tran- 
sitions due to the large linewidths observed toward the GC 



clouds. To account for the overlap, we used the following con- 
strains for the Gaussian fitting depending on the observed tran- 
sitions: 

- For molecules without internal rotation, like CS, C ls O, 
13 CO, HCOOH, and H' 3 CO, and for molecules with internal 
rotation when the transitions are sufficiently separated, like 
those of C2H5OH and some of the HCOOCH3 lines, simple 
Gaussian profiles were fitted to each transition for every 
velocity component. As mentioned before, the CS, C 18 0, 
and 13 CO profiles generally show more velocity components 
than those observed in the complex organic molecules. For 
this work we only fitted the components in the velocity range 
where at least one organic molecule was detected. In the case 
where the velocity components overlapped we calculated the 
integrated intensity of CS, C I8 0, and 13 CO in the velocity 
range where the complex organic molecules show the emission. 

- The CH3OH and I3 CH30H profiles are composed of 
several overlapped transitions arising from levels at different 
energies. In this case, we fitted the blended transitions with 
multi-Gaussian profiles forced to have the same linewidth and 
velocity separations corresponding to the rest frequencies of 
the transitions. We also forced the relative intensities derived 
from the spectroscopic parameters (Tables [0 and [3}, assuming 
the same rotational temperature for all the transitions. In 
addition to the radial velocity and linewidth, we also fitted 
the rotational temperature and the optical depths that match 
the observed profiles. In this case, the antenna temperatures 
were obtained from (Trot-TbgXl-e -1 "). We used the optically 
thin emission from I3 CH30H transitions to test the optical 
depths and the antenna temperatures obtained from the 
CH3OH fits. The optical depths derived from the CH3OH 
lines are in good agreement with those obtained from the 
[CH3QH/ 13 CH30H1 line-intensity ratios, assuming a 12 C/ 13 C 
ratio of 20 JWilson & Roodlll994 . For the 13 CH 3 OH, and in 
some cases for the CH3OH lines, only the transitions with 
the largest intensities were fitted. In the case of the CH3OH 
5o — > 4o£' transition, a simple Gaussian fit was made since this 
line is not overlapped with any other transition. 

- Due to its internal rotation, the (CITO2O rotational levels 
are split into four substates with similar Einstein coefficients 
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Table 2. Source position 



Requena-Torres et al.: Organic molecules in the GC 



source 


1(°) 


bC) 


a (B 1950) 


6 (B1950) 


Region 


MC G-0.96+0.13 


-0.96 


+0.13 


1? h 39m 


36 s 6 


-29° 39' 47" 


SgrC 


MC G-0.55-0.05 


-0.55 


-0.05 


17 h 


41 m 


20 s 


-29° 24' 30" 


SgrC 


MC G-0.50-0.03 


-0.50 


-0.03 


17 h 


41 ra 


21 s 2 


-29° 21' 26" 


SgrC 


MC G-0.42-0.01 


-0.42 


-0.01 


17 h 


41 ra 


24 s 1 


-29° 15' 51" 


SgrC 


MC G-0.32-0.19 


-0.32 


-0.19 


17 h 


42 m 


24 s 7 


-29° 17' 20" 


SgrC 


MC G-0. 11-0.08 


-0.11 


-0.08 


17 h 


42 m 


28 s 


-29° 02' 55" 


SgrA 


MC G-0.08-0.06 


-0.08 


-0.06 


17 h 


42 m 


30 s 


-29° 00' 58" 


SgrA 


MC G-0.02-0.07 


-0.02 


-0.07 


17 h 


42 m 


40 s 


-28° 58' 00" 


SgrA 


MC G+0.04+0.03 


+0.04 


+0.03 


17 h 


42 m 


26 s 2 


-28° 51' 45" 


TRA 


MC G+0.07-0.07 


+0.07 


-0.07 


17 h 


42 m 


54 s 2 


-28° 53' 30" 


TRA 


MCG+0. 13+0.02 


+0.13 


+0.02 


17 h 


42 m 


41 s 4 


-28° 47' 35" 


TRA 


MCG+0. 17+0.01 


+0.17 


+0.01 


17 h 


42 m 


50 s 


-28° 45' 50" 


TRA 


MCG+0. 18-0.04 


+0.18 


-0.04 


17 h 


43 ra 


0P0 


-28° 47' 15" 


Sickle 


MC G+0.20-0.03 


+0.20 


-0.03 


17 h 


43 m 


03 s 6 


-28° 45' 42" 


Sickle 


MC G+0.24+0.01 


+0.24 


+0.01 


17 h 


42 m 


59 s 6 


-28° 42' 35" 


Sickle 


MCG+0.62-0.10 


+0.62 


-0.10 


17 h 


44 m 


18 s 


-28° 26' 30" 


Sgr B2_a 


MC G+0.64-0.08 


+0.64 


-0.08 


17 h 


44 m 


17 s 5 


-28° 24' 30" 


SgrB2_a 


MC G+0.67-0.06 


+0.67 


-0.06 


17 h 


44 ra 


18 s 


-28° 22' 30" 


SgrB2_a 


MC G+0.68-0.10 


+0.68 


-0.10 


17 h 


44 ra 


27 s 2 


-28° 23' 20" 


SgrB2_a 


MC G+0.70-0.01 


+0.70 


-0.01 


17 h 


44 ra 


10 s 


-28° 19' 30" 


Sgr B2_a 


MC G+0.70-0.09 


+0.70 


-0.09 


17 h 


44 ra 


27 s 2 


-28° 22' 05" 


Sgr B2_a 


MCG+0.71-0.13 


+0.71 


-0.13 


17 h 


44 ra 


38 s 4 


-28° 22' 25" 


Sgr B2_a 


MC G+0.76-0.05 


+0.76 


-0.05 


17 h 


44 ra 


27 s 2 


-28° 17' 35" 


Sgr B2_a 


SGR B2N 


+0.68 


-0.03 


17 h 


44 ra 


10 s 6 


-28° 21' 17" 


Hot Core 


SGR B2M 


+0.66 


-0.04 


17 h 


44 ra 


10 s 6 


-28° 22' 05" 


Hot Core 


MC G+0.694-0.017 


+0.694 


-0.017 


17 h 


44 m 


10 s 


-28° 20' 05" 


SgrB2 


MC G+0.693-0.027 


+0.693 


-0.027 


17 h 


44 ra 


12 s 1 


-28° 20' 25" 


SgrB2 


MC G+0.627-0.067 


+0.627 


-0.067 


17 h 


44 m 


12 s 1 


-28° 25' 05" 


SgrB2 


MC G+0.630-0.072 


+0.630 


-0.072 


17 h 


44 m 


13 s 6 


-28° 25' 05" 


SgrB2 


MC G+0.672-0.014 


+0.672 


-0.014 


17 h 


44 ra 


06 s 1 


-28° 21' 05" 


SgrB2 


MC G+0.640-0.046 


+0.640 


-0.046 


17 h 


44 m 


09 s 1 


-28° 23' 45" 


SgrB2 


MC G+0.635-0.069 


+0.635 


-0.069 


17 h 


44 m 


13 s 6 


-28° 24' 45" 


SgrB2 


MC G+0.659-0.035 


+0.659 


-0.035 


17 h 


44 m 


09 s 1 


-28° 22' 25" 


SgrB2 


MCG+0.68 1-0.028 


+0.681 


-0.028 


17 h 


44 m 


10 s 6 


-28° 21' 05" 


SgrB2 


MC G+0.673-0.025 


+0.673 


-0.025 


17 h 


44 ra 


09 s 1 


-28° 21' 25" 


SgrB2 



Note. Galactic and equatorial coordinates of the selected sources. Nomenclature, as recommended by the IAU: MC followed by a G and the 
galactic coordinates, where MC stands for Galactic center molecular cloud. We have grouped the sources by regions: sources around Sgr C, 
Sgr A, the TRA, the Sickle, Sgr B2 (_a is to difference two observation sessions), and the hot cores. Not all sources in the TRA or in the Sickle 
are affected by the same physical conditions. 



but different intensities due to the nuclear spin degeneracy. 
These substates are called AA, EE, EA, and AE. Since the four 
substates are blended, we measured the integrated intensity 
of the four substates and then summed all the spin weights of 
each substate to derive the (CH3)20 column densities. 

- To obtain the upper limits to the column densities of 
undetected transitions like those of CH3COOH, we used the 
3<x level for the integrated intensities. 

Derived parameters from the Gaussian fits (peak intensity, 
integrated intensity, radial velocity and linewidth) for all transi- 
tions in each source are shown electronically in Table ??, avail- 
able at the CDS. 

The linewidths of the complex organic molecules in our 
sources are ~ 15 km s _1 , typical in the GC region. Radial veloc- 



ities of the different sources range from -93 km s to 140 km 

s" 1 . 



4. Derived parameters 

4.1. Column densities, excitation temperatures, and 
densities 

Molecular column densities were derived by assuming opti- 
cally thin emission and the local thermodynamic equilibrium 
(LTE) approximation. In the case of optically thick lines like 
those of CH3OH, the optically thin isotopic substitution was 
used. Under these conditions, the total column density, N, of a 
molecule is given by: 

N = 3L Qe (EJ«T M ) (1) 
gugi 



Table 3. Dipole moments and rotational constants of the observed molecules 



molecule 






A 


B 


C 


Q 


comments 


references 




(Debye) 


(Debye) 


(MHz) 


(MHz) 


(MHz) 








CH3OH 


0.885 


1.440 


127484.00 


24679.98 


23769.70 


1 28 r 1 ? 


Rotational levels split into E and A substates. 
Observed transitions with several blended lines. 


Menten et al. (1986) 
Anderson et al. (1990a) 


13 CH 3 OH 


0.899 


1.440 


127527.40 


24083.50 


23197.00 


1.28 Tg 


Rotational levels split into E and A substates. 
Observed transitions with several blended lines. 


Anderson et al. ( 1990b) 
Xu & Lovas (1997) 


C2H5OH 


0.046 


1.438 


34891.77 


9350.68 


8135.34 


3.2 rif 


Rotational levels split into Trans and Gauche substates. 
Only Trans Transitions are exited at ~ 10 A". 


Pearson et al. (1995. 1997) 


HCOOCH3 


1.630 


0.680 


19985.71 


6914.75 


5304.48 


12.45 


Rotational levels split into E and A substates. 


Plummer et al. (1984. 1986) 
Kuan & Snvder (1996) 


(CH 3 ) 2 


0.000 


1.302 


38788.20 


10056.50 


8886.80 


91.6 T£ 


Rotational levels spilt into AA, EE, AE and EA substates. 
g,: (K n K r : ee <-> 00) AA:EE:AE:EA=6: 16:2:4 
gr. (KaK^: oe <-» eo) AA:EE:AE:EA= 10: 16:6:4 


Groneretal. (1998) 


HCOOH 


1.396 


0.260 


77512.25 


12055.11 


10416.12 


1.68 Tg 




Willemot et al. (1980) 














14.1 r r ' D f 




Liu et al. (2002) 


CH3COOH 


0.860 


1.470 


11335.58 


9478.73 


5324.99 
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Fig. 2. Observed line profiles toward the sources MC G-0.02-0.07, MC G+0.76-0.05, and MC G+0.20-0.03. Line intensi- 
ties are expressed in antenna temperature (K), and velocities refer to LSR (km s _1 ). The CH3OH profiles are composed of 
several overlapped transitions arising from levels with different energies; in the MC G-0.02-0.07 source, we show the tran- 
sitions with highest intensities as vertical solid lines and the transitions with lowest intensities as dotted lines. These transi- 
tions are: forCH 3 OH J=3->2, from left to right, 3 2 ^2 2 A, 3i^2iE, 3- 2 ^2_ 2 E, 3 2 ^2 2 E, 3 2 ^2 2 A-, 3 ^2 A+, 3-i^2_iE, and 
3 ^2 E. For the 13 CH 3 OH J=3^2, from left to right, 3 2 i->2 20 A+, 3i 2 ->2 n E, 3 2 i->2 21 E, 3_ 22 ->2_ 21 E, 3 22 ^2 2i A-, 3 03 ->2 02 A+, 
3_i 3 ->2_ 12 E, and 3 03 ->2 02 E. For the CH3OH J=5->4 from left to right; 5 2 ->4 2 E, 5_ 2 ->4_ 2 E, 5 2 ^4 2 A+, 5i^4iE, 5_ 3 ->4_ 3 E, 
5 3 ->4 3 E, 5 2 ->4 2 A-, 5 3 ->4 3 A-, 5 3 ^4 3 A+, 5 4 ^4 4 E, 5_ 4 ^4_ 4 E, 5 4 ^4 4 A-, 5 4 ->4 4 A, 5 ^4 E, 5^ ->4_iE, and 5 ->4 E. All these 
transitions have been used for the multi-Gaussian fit. The spectrum toward MC G+0.20-0.03 shows two velocity components. 



N u = 1.67 xlO 14 ^- T' A dv, (2) 
vu l S J 



where N„ is the column density in the upper level of the 
observed transition in crrT 2 , g u the upper level degeneracy 
(2J U +1), v the line frequency in GHz, /u the dipole moment in 
Debye, S the line strength, J T A dv the integrated intensity in K 
km s -1 , gi the spin weight degeneracy, Q the partition function, 
E u /k the upper level energy in K, and T rot the rotational tem- 
perature also in K. Table^summarizes the dipole moments, the 
rotational constants, the partition functions, and the statistical 
weights used to calculate the total column densities. 



4.1 .1 . Rotational temperatures 

To derive the total column densities one first needs to esti- 
mate T Iot . For molecules with more than one observed tran- 
sition, we calculated T I0t from population diagrams. Figure^] 
shows, as an example, the population diagrams (ln(N u /g u ) 
vs. E u /k) of some of the organic complex molecules derived 
for the source MC G-0.02-0.07. In general, r rot from differ- 
ent molecules (third column on Table [6} agree within a fac- 
tor of 2. We find r, ot values between 6K and 18K, much 
lower than the typical kinetic tempera t ures of ^:100K de- 
rived from NH 3 bv iHutte meister ■ e\ al.| dl993al) and H 2 by 
iRodriguez-Fernandez et al.l J20001 l2001al) . indicating subther- 
mal excitation. We do not find any systematic trend in r rot as a 
function of the molecule. For sources with only one observed 
transition, we assumed an average value of T rot ~8 K. 
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Fig. 3. Same as Fig. [2] The 13 CO, C 18 0, and CS profiles show wider lines than those of the complex organic molecules in Fig. [2] 



4.1 .2. Densities 

Once the T mt is derived, we can make a rough estimate of the 
H2 densities, providing that the kinetic temperatures and the 
collisional cross section are known. The H2 densities can be 
estimated from 



1 + V ~ ex P{7t)\ 



(3) 



where A„i is the Einstein coefficient between the upper and 
the lower levels, <v<x> u i the collisional rate coefficient be- 
tween the same levels, v the transition frequency, h the Plank 
constant, k the Boltzmann constant, T mt the rotational tem- 
perature, Th„ the background tem perature, and 7\ the kinetic 
temperature. Pottag e et alJ ll2004 have calculated the colli- 
sion rate coefficients (<vcr>) between E-CH3OH and para-H2 
for kinetic temperatures between 5 K and 200 K. As suming 
a kinetic temperature of 100 K fHtittemeist er et all [l993a; 

nguez- Fernandez et"ai1l2001ah . we obtain H2 densities of 
~10 5 cm 3 for the rotational temperature derived from CH3OH 
of 6-18 K. Similar H2 densities are also derived when the mea- 
sured line intensities of E-CH3OH are fitted using the large 
velocity gradient (LVG) approximation model, provided by J. 
Cernicharo, for the excitation of E-CH3OH. We also estimated 
the H2 densities from the excitation temperatures derived from 
C2H5OH. The collisional cross section in this case was esti- 
mated from those of CH3OH but corrected for the difference in 



mass and size. These roughly estimated H2 densities are also 
~10 5 crrT 3 . 

The derived H2 densities are based on the collisional rates 
for a kinetic temperature of 100 K. However, the linewidths in 
the GC molecular clouds indicate the presence of supersonic 
turbulence and thus of non-equilibrium processes. In this case 
the collisional rates will not be described by the kinetic tem- 
peratures. However, for the observed transitions the rate coef- 
ficie nts do depend w eakly on the increasing collisional veloci- 
ties JPottage et all2004l) . therefore our estimated densities will 
only be marginally affected by the effect of a larger velocity dif- 
ference between the collisional partners than those represented 
by the kinetic temperature. 

4.1 .3. The H 2 column densities 

To estimate the H2 column density, we used the J=l— >0 and 
J=2^1 C 18 and the J=l->0 13 CO transitions using the LTE 
approximation . We used the isotopic ratios of 12 C/ 13 C= 20, 
16 0/ 18 0= 250 JWilson & Roodll994l). and a CO relative abun- 
dance of CO/Ht = 10~ 4 jFrerking et alJl982t) . The last column 
in Table[6]shows the derived H2 column densities. Our H2 col- 
umn densities, derived from the LTE approximat ion, are w ithin 
a factor of <2 of those derived by Rodriguez-Ferna ndez et alJ 
ll2001ah using the LVG approximation. 

Table [6] also shows the relative CS abundance estimated 
from the CS column densities derived using the T mt obtained 
from CH3OH. Since our data refer to the line emission of the 
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Fig. 4. Population diagrams for some complex molecules to- 
ward MC G-0.02-0.07. The error bars are derived by assum- 
ing a calibration uncertainty of 10% and the Gaussian fit error. 
Parentheses () indicate the calculated errors. Triangles with ar- 
rows represent upper limits to HCOOCH3 column densities. 



main CS isotopomer, optical depth eff ects migh t affect the de- 
rived CS column densities dHuttemeisterll993b|). We compared 
our results with those derived from 13 CS by Mart in et al 1 d2006l) 
in several common sources. The CS column densities derived 
from I3 CS are always higher than those of the CS by a factor 
of 1-5. When the 13 CS determinations are available, Table |6] 
shows the two values estimated for the relative abundance of 
CS. 



4.1 .4. Column densities of complex organic molecules 

To derive the column densities of the complex organic 
molecules we used the parameters in Table [5] and Eqs.^and 
|3 The relative abundance of complex molecules derived from 
their column densities and the H2 column densities are shown 
in Table |6] 



4.2. Abundances and abundance ratios 

Table |6] shows the derived rotation temperatures, fractional 
abundances, and the H2 column densities for the different 
molecular clouds. Different sources in the same line of sight 
are identified by their radial velocities. The fractional abun- 
dances of the complex organic molecules are very high, even 
higher than in the hot cores. The abundances we find reach val- 
ues up to 1.1x10" forCH 3 OH, 6xl0~ 8 forC 2 H 5 OH, 5.6xl0~ 8 
for (CH 3 ) 2 0, 7.5 xl0~ 8 for HCOOCH3, 4.1X10 -9 for HCOOH, 
and 1.9xl0~ 8 for H2CO. We have only obtained upper lim- 
its to the abundance of CH3COOH between <5xl0~ 10 and 
<5.5xl0~ 8 . 

There are some sources for which only upper limits to 
the abundances were measured for all the complex organic 
molecules except for CH3OH. In some cases, like the sources 
with galactic longitude <-0.30°, the upper limits to the abun- 
dances of all complex molecules are not relevant because of 
the lack of sensitivity of our observations. In this case the 
upper limits have not been included in the discussion. Other 
cases, like the C 2 H 5 OH upper limits toward MC G+0. 18-0.04, 
MC G+0.20-0.03 (Sickle), and MC G+0. 13 +0.02 (TRA) are 
relevant and therefore included in our discussion. The molecu- 
lar emission toward the Sickle shows two velocity components 
at ~25 km s _1 and ~80 km s _1 , which are believed to be af- 
fected by the UV radiation p roduced by the Quintuplet cluster 
dRodriguez-Fernandez et al]l2001bl) . The molecular emission 
toward MC G+0. 13 +0.02, located in the TRA, also seems to 
be affected by the UV radiation from the Arches Cluster. In fact 
only the -30 km s _1 velocity component, which is clear ly asso- 
ciated with the ionized gas dSerabvn & Giisten1ll987l) . shows 
lower CaHgOH abundance than in other sources. 

Because CH3OH is the only molecule that has been 
detected in all sources and it is believed to play a central role 
in the formation of some of the observed complex organic 
molecules, we compared the abundances of the rest of the 
complex molecules with that of CH3OH. Figure|5]summarizes 
the results illustrating the abundance of all the observed 
molecules relative to that of CH3OH as a function of the 
CH3OH abundance. As already mentioned, we have only 
included the significant upper limits in this figure. The main 
results are: 

-In spite of the C^OH's abundance changing by a factor 
of -50, all the complex organic molecules, except that of the 
CH3COOH with only upper limits to its abundance, show 
a surprisingly constant relative abundance with respect to 
CH3OH. 



Requena-Torres et al.: Organic molecules in the GC 



9 



'0 



o 

X 

P 10 



P 10 



GC 



i SgrB2 
SgrB2 a 
SgrA 
Sickle 
TRA 



. «:. v 



! t 



t 1 



Tr 



CHjOH/H 2 



Fig. 5. Relative abundances of organic complex molecules with 
respect to CH3OH as a function of the CH3OH abundance 
for the selected GC clouds. Different regions are represented 
with different symbols (see Table|2ji. The arrows represent up- 
per limits. Only the detections and significant upper limits are 
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CH 3 0H/H 2 

Fig. 6. Same as Fig.|5]but for H 2 CO. 

-The mean [C2H5OH/CH3OH] abundance ratio of 
~3.6xl0~ 2 is constant within a factor 3.7. The Sickle and 
the TRA clouds show upper limits to the [C 2 H 5 OH/CH 3 OH] 
ratios, which are a factor of k,6 times lower than the mean ratio 
found in the typical GC clouds, a factor of SI. 6 higher than 



the scattering of the data. 

-The abundance ratios [(CH 3 ) 2 0/CH 3 OH] and 
[HCOOCH3/CH3OH] present similar mean values of 
~3.9xl0~ 2 . These ratios show a somewhat larger disper- 
sion than that of C2H5OH but are constant within a factor of 7. 

-For the [HCOOH/CH3OH] ratio, we obtained a mean 
value of ~7.9xl(T 3 and a dispersion of a factor of 5.3. 

-CH3COOH was not been detected in the GC clouds but 
its [CH3COOH/CH3OH] ratio must be lower than 10 3 . 

-We also show in Fig. |6] the [H 2 CO/CH 3 OH] ratio as a 
function of the CH3OH abundance. As for the more com- 
plex organic molecules, we find a constant ratio of ~l.lxl0~ 2 , 
within a factor 3. 
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Fig. 7. Variation of [CS/CH3OH] as a function of CH3OH 
abundances. The vertical bars show the changes due to the 
opacity. 

The [CS/CH3OH] abundance ratio is shown in Fig.0 This 
abundance ratio decreases by a factor 60 as the CH3OH abun- 
dance increases. This is the only molecule for which we find a 
systematic trend, and the variation is clearly larger than those 
found for the complex organic molecules. As previously men- 
tioned, our CS column densities could be affected by opac- 
ity effects. To show the effects of opacity, we have also in- 
cluded the abundance ratio obtained using the rarer CS iso- 
topomer 13 CS in Fig. |7] as vertical bars. The [CS/CH3OH] 
abundance ratio derived from the optically thin isotopomer in- 
creases by a relatively constant value for all clouds indepen- 
dently of the CH3OH abundance. The variation observed in the 
[CS/CH3OH] ratio as a function of the CH3OH abundance, us- 
ing the 13 CS data, is similar or even larger than that derived 
from the CS data. We therefore consider that the difference be- 
tween the behavior of CS and the complex organic molecules 
is real and related to the different type of chemistry in the GC 
clouds. 



5. Comparison with Galactic disk objects 

In dark clouds, where the kinetic temperature is only 10K, the 
abundance of complex organic molecules is very low. In the 
dark cloud TMC-1, the CH3OH, H 2 CO, and HCOOH abun- 
dances are as low as 2x1 (T 9 , 2xl(T 8 , and 3xl0 -10 , respec- 
tively (see e.g. lTurnerl2000l) . For other dark clouds like L183, 
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one finds similar abundances of 8x10 9 for CH3OH, of 2x10 8 
for H,CO, and of 3xl0 10 for HCOOH dTurner et al.lll999t 
iDickens et alT2 000). Other complex molecules, like CaHgOH, 
(CH 3 ) 2 0, CH3COOH, and HCOOCH3, have not been detected 
so far, but model pr edictions suggest C2H5OH abundances of 
10" 12 (flasegawa & Herbst fl993l) . High abundance of the com- 
plex organic molecules were first observed in the short-lived 
(10 5-6 years) objects associated with massive star-forming re- 
gions known as hot cores. There, it is believed that grain 
mantle evaporation relea se complex mol e cules, like alco hols, 
into gas phase. R ecently, iBottinelli etaD d2004albl l200fj) and 
lKuanetalJll2004) detected relatively large abundances of com- 
plex organic molecules, like HCOOCH3 and (CH^O, toward 
the so-called "hot corinos" associated with warm condensation 
surrounding low mass proto-stars. 

In the following sections, we compare the abundances of 
complex organic molecules in the GC with those found in the 
hot cores and hot corinos associated with star formation. 



5.1. Hot cores 

For the comparison between the GC and hot cores associated to 
mass ive star formatio n, we selected the homogeneous data set 
bv llkeda et all d20Qlh for several hot cores (W51, NGC6334f, 
G327.3-0.6, G31.41+0.31, G34.3+0.2, G10.47+0.03, and the 
Orion Hot Core). We also used data from lRemiian et alJ (12002. 
120031) to obtain the CH3COOH abundances derived from their 
[HCOOCH3/CH3COOH] ratio. The hot cores abundances of 
complex molecules and the H2 column densities are given in 
Table 0] Both the H2 and molecular column densities were 
derived from observations with similar beam sizes. Then the 
derived fractional abundances are independent of beam dilu- 
tion. Since our comparation is based on the relative abundance 
between complex molecules, our results are not affected by the 
size of the hot cores. We also used our data for the Sgr B2N 
and Sgr B2M hot cores in Table [6] Figure [8] illustrates the 
comparison between the abundance ratios observed for the 
GC and those measured in hot cores. As for the GC clouds, 
the abundance of CH3OH in hot cores shows variations of 
more than one order of magnitude between different objects. 
Surprisingly, the abundance ratios relative to CH3OH are also 
independent of the CH3OH abundance and similar to those 
measured for the GC clouds. The C2H5OH and HCOOH mean 
abundance ratios with respect to CH3OH between hot cores 
and GC clouds are different by less than a factor of 2. The 
molecule that presents the largest difference in its abundance 
ratio with respect to that in the GC clouds is (CH3)20 with a 
mean ratio of 2xl0~', which is ~6 times higher than those of 
the GC clouds. The [CH3COOH/CH3OH] ratios in hot cores 
agree with the upper-limit ratios found in GC clouds. Our data 
for the Sgr B2N Hot Core shows slightly higher abundances 
of HCOOCH3 and HCOOH than the rest of the hot cores. The 
Sgr B2M Hot Core, with the lowest CH3OH abundance, shows 
similar ratios to those in the GC sources. 

It is remarkable that different objects like the hot cores and 
the GC clouds, with different properties in terms of densities, 
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Fig. 8. Same as Fig. |5 ] bu t inclu d ing the hot core abundance s 
obtained from llkeda et alJ d200ll) . lRemiian et all J2002il2003h . 
and from our SGR B2 data as open symbols. 



H2 column densities, and spatial extension, show such uniform 
abundance ratios of very complex molecules. Hot cores are 
denser than the GC clouds by more than 2 orders of magni- 
tude, and they typically have larger column densities than the 
GC clouds. 



5.2. Hot corinos 



IBottinelli et all d2004albl l200d) . Lforgensen etaD d2005l) 
and Ku an et alJ ( 120041) observed the hot corinos associ- 
ated with the low mass proto-stars, NGC1333 IRAS4A, 
NGC1333 IRAS4B, NGC1333 IRAS2A and IRAS 16293- 
2422 (A and B), which also show large abundances of 
HCOOCH3, H 2 CO, (CH 3 ) 2 and HCOOH. Table g] summa- 
rizes the recent compilation o f the hot corino abundances of 
complex organic molecules bv IBottinelli et alJ d2006l) . In this 
case the molecular abundances were derived by as suming a 
source size of 0.5" for IRAS 4A dBottinelli et al l2004al) . 0.25" 
for IRAS 4B, 0.43" fo r IRAS2A dMaret et alJl2004 . and 2" 
for IRAS 16293-2422 dCazaux et al.l2003h . 
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Except for IRAS 4A, which shows lower abundances, 
HCOOCH3, (CH 3 ) 2 0, HCOOH, and H 2 CO show higher abun- 
dances in hot corinos than in the GC molecular clouds and hot 
cores. The data show that the abundance of CH3OH varies be- 
tween sources by nearly two orders of magnitude, similar to the 
variation observed in the GC clouds. As in the GC molecular 
clouds and in the hot cores, the abundance of the complex or- 
ganic molecules relative to that of CH3OH is independent of the 
CH3O H abundance for different hot corinos (Bot tinelli et alJ 
2006). 

Although the number of observations of hot corinos is still 
limited, it seems that the abundance ratios, including species 
believed to be formed by different chemistries, show a sim- 
ilar behavior to the GC sources. However, the abundance of 
the molecules is very different, as expected if the grain man- 
tle composition in dark clouds were different than in the GC 
clouds and hot cores. 



6. Discussion 

With the exception of the clouds associated with the Sickle 
and the TRA, large abundances of organic complex molecules 
like CH3OH, C 2 H 5 OH, (CH 3 ) 2 0, HCOOCH3, HCOOH, and 
H 2 CO seems to be a general characteristic of molecular clouds 
in the GC at scales of a few hundred parsecs. The main 
characteristic of the GC molecular gas is that the complex 
organic molecules show a constant relative abundance with 
respect to CH3OH, within a factor 4-8 over the whole range 
of the CH3OH abundances in the observed GC clouds. This 
contrasts with the [CS/CH3OH] abundance ratio that increases 
by a factor of 60 when the CH3OH abundance decreases. 

We now discuss the origin of the chemistry of the com- 
plex molecules found in the GC sources, including the low 
[C 2 H50H/CH30H] abundance ratio measured in the Sickle and 
the TRA clouds. 



6.1. Models for gas-phase formations of organic 
molecules 

Table|5]shows a compilation of the models proposed to explain 
the formation of complex organic molecules. Since gas-phase 
chemistry cannot account for the large abundance of alcohols, 
such as CH3OH and C 2 H50H, or for other organic molecules 
like the H 2 CO measured in hot cores, evaporation from grain 
mantles has been proposed as the main formation mechanism. 
Most of these models assume that (CH 3 ) 2 and HCOOCH3 
were produced in gas phase after the ejection of CH3OH. 
However, it has been recently claimed that HCOOC H3 cannot 
be pro duced efficiently in gas phase from CH3OH dHorn et alJ 
2004). The large abundances of this molecule found in the GC 
clouds clearly indicate that HCOOCH3 is hardly likely to be 
produced in gas pha se. However, the most recent model from 
IPeeters et alJ J2006I) studied the possibilities for the formation 
of (CH3) 2 on grain mantles and concluded that the gas-phase 
reactions are the dominant way to form (CH3) 2 0, while grain 
mantle reactions are the minor source. 

From the theoretical point of view, so far, only (CH3) 2 
seems to be produced in gas phase from reactions involving 
CH3OH on relatively short time scales. Assuming that the 
gas phase models for the (CH3) 2 formation are correct 
(see references in Table one would expect large changes 
in the relative abundances between the daughter ((CH3) 2 0) 
and parent (CH3OH) molecules as a function of time after 
ejection/evaporation. One expects low [(CH3) 2 0/CH30H] 
abundance ratio in early times just when large abundances of 
CH3OH are ejected into the gas phase and the gas phase pro- 
cessing starts. Low [(CH3) 2 0/CH30H] ratios are also expected 
in late times in the low CH3OH abundance regime when an 
important fraction of CH3OH has been converted into daughter 
molecules and both molecules have been destroyed by gas- 
phase reactions. Between these two regimes, one would expect 
the [(CH3) 2 0/CH30H] ratio to reach its maximum values, be- 
cause the daughter molecules reach their maximum abundance. 
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Table 5. Summary of chemical models and observations of complex organic molecules 



source 


T 

(K) 


n(H 2 ) 

(cm- 3 ) 


CiH 5 OH n 
CHiOH 

xlO" 2 


CH,0H 

xlO" 2 


HCOOL'th a 
CHiOH 

xlO" 2 


Hi CO a 
CH 3 OH 


age 
xl0 5 y 


Galactic center clouds 


-100* 


~10 4 -10 5 


-3.66 


-3.87 


-3.96 


-LlOxlO 2 




hot cores r 


-100 


~10 6 


-6.90 


-21.80 


-8.03 






Leung et al. (1984) rf 


10 


10 4 


0.14 


6.00 




43.6 


3.16 


Herbst & Leung (1986) e 


10 


10 4 


2.18 


5.25 




9.75 


3.20 


Herbst & Leung (1989)'' 


10 


10 4 


0.99 


0.12 




65.5 


1.00 


Millar et al. (1991V 


70 


2xl0 5 


0.19 


6.32xl0- 2 


1.26 


1.18xl0 3 


1.60 


Hasegawa & Herbst (1993/ 


10 


10 4 


0.77 


1.56xl0 2 




4.06xl0 2 


3.20 


Caselli et al. (1993)-CR S 


100 


10 6 


1.5X10" 4 


0.21 


4.67x10-3 


1.05 


1.00 


Charnlev et al. (19951 


100 


2xl0 6 




-2 


-1.5 


4x1 0~ 2 


-0.30 


Rodgers & Charnlev (2001)'' 


100 


10 7 


3.75 


30.00 


10.00 




0.12 


Rodgers & Charnlev (2001V 


100 


10 7 


20.00 


12.50 


0.50 




0.60 


Horn et al. (2004V 


100 


10 6 




0.75 


3.33xl0- 2 


1.56 


0.80 


Peeters et al. (2006)* 


100 


10 7 




50.00 






0.20 



" The relative abundances of C 2 H 5 OH, (CH 3 ) 2 0, HCOOCH3, and H 2 CO with respect to CH3OH predicted by time-dependent chemical 
models, r efer to the age of m aximum abun dance of the CH 3 OH's dau ghters. 
b From Hi ittemeister et alj Jl993al) and Rodriguez-Fernand ez et alJ ^2001a). 
Mean abundance ratios for the hot core sources used in this work. 
d Model 1 of the referred paper. 
e Low-metal model in the referred paper. 
f Model starting with all the H in H 2 . 
8 Model fixed with the Orion Compact ridge observations. 
* Model without NH 3 . 
' Model with NH 3 . 
' Model 3 of the referred paper. 
k Model with (CH 3 ) 2 created only in gas phase. 



From the hot core models, one expects variations in 
[(CH3)20/CH30H] abundance ratio up to three orders of 
magnitude as a function of time. This contrasts with our data 
for the GC clouds and the hot cores, which show variations of 
less than one order of magnitude. Assuming that the chemical 
models for hot cores are applicable in the GC, the only possi- 
bility of explaining the high abundances of complex molecules 
and the rather constant [(CHs^O/CH^OH] abundance ratios 
is to consider that all GC clouds in our study, distributed 
over 200 pc, have undergone the ejection of CH3OH from 
grain mantles nearly simultaneously, -10 5 years ago. This 
possibility seems very unlikely, suggesting that (CH3)20 is 
also ejected from grain mantels. Since (CtT^O shows similar 
abundances with respect to CH3OH than the other molecules 
that can only be formed on grain mantles like C2H5OH and 
HCOOCH3, our data support the scenario in which all complex 
organic molecules in the GC clouds, the hot cores, and also in 
hot corinos have been ejected/evaporated from grains. 



Previously we also found differences in the behavior of the 
relative abundances of CS and the complex organic molecules. 
This behavior is supported by t he data in shock regions as soci- 
ated with outfl ows, as in LI 44 8 ( Jimenez-S erra et al.l2 005) and 
LI 157 (Bachi ller & Perez Gutierrez! 19971) . The CS abundance 
is enhanced by the shock, but the enhancement of CH3OH is ~ 1 
order of magnitude larger than that of CS. This effect of differ- 
ential ejection of the molecules can explain the variation in the 



complex organic molecules abundances in gas phase, while the 
CS abundance remains marginally affected. 

6.2. Grain mantle erosion by shocks and time scales 
for depletion 

In the GC clo uds, the typical dust-grain tempe ratures are low 
(T~ 10-30 K, Rodngue z^Fernandez et"ai1 E004) and the ejec- 
tion of molecules to gas phase cannot be produced by evap- 
oration like in hot cores, but by shocks with moderate ve- 
locities. The shock heats the gas and produces the sputter- 
ing of molecules on the grain mantles. For shock veloci- 
ties of ~20km s" 1 , the ga s temperature can reach -1000K 
dKau fman & Neufeldlll996l) . After the shock passage, the gas 
rapidly cools on time scales of ~10 4 years for densities of 
10 5 cm 4 . The complex molecules ejected from the grains will 
again stick to the grains or grain mantles. The time scales 
will depend on the sticking coefficient. A sticking coefficient 
of ~1 is expected for mean velocities corresponding to the 
thermal velocity at temperatures of <200 K because the com- 
plex molecules will be absor bed after every collision with a 
grain dTielens & Hagenlll982l) . However, the sticking coeffi- 
cient decreases when increa sing the relative velocity betwee n 
the grains and the mol ecul es jLeitch- Devlin & WilliamslT 9 85). 
Buch & Zhand dl99ll) and lMasuda et alJ dl998l) estimatedthat 
the sticking coefficient of H is a function of the mean velocity. 

The key parameter is the adsorption energy of the com- 
plex molecules, which will depend on the surface of the dust 
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grains. Since the CH3OH abundance in gas phase is rela- 
tively large, grain mantles are completely eroded and the most 
likely grain surfaces in the GC clouds would be graphite and 
silicates. The abso rption energy of CH3 OH on carbon basal 
planes is 1600 K (Aikawaetal. 1996) and on silicates is 
2065 K ilAllenfcRohinsorJll977h . For the typical mean ve- 
locities in the GC of 2-4 km s (see next section) the stick- 
ing coefficient will decrease and the time scale for sticking 
would be ~10 5 years on silicates and graphites. Similar time 
scales are found for the depletion of C2H5OH, HCOOCH3, and 
(CH0? O molecules with absorp tion energies ~3000 K on sil- 
icates dAllen & Robinsonlll977l) . Short time scales for deple- 
tion, like those expected in the GC, would require a mechanism 
that continuously ejects molecules from the grain mantles. 

6.3. Origin of the shocks in the GC 

Rodriguez-Fernand ez et al ] (bom id )2004) have summarized 
the possible mechanisms for the heating and the chemistry 
of the molecular gas in the GC. Turbulence and low-velocity 
shocks are proposed as the most likely mechanism for driving 
the chemistry. One of the key parameters for constraining the 
origin of the shocks is the velocity dispersion in the GC molec- 
ular clouds. High angular-resolution observation of the enve- 
lope of Sgr B2 indicates that the warm molecul ar gas is hi ghly 
turbul ent with line widths of ~4km s" 1 jMartin-Pin tado et alJ 
h999USrrrithetal]lb000h show that MHD turbulence creates a 
wide range of shock velocities, but the larger amount of them 
should be produced with Mach numbers between 2-4. In the 
GC regions the sound speed is ~1 km s _1 (at 100 K), and then 
the more abundant shocks associated with turbulence should 
have velocities between 2 and 4 km s -1 , in agreement with the 
linewidth measured from high angular-resolution observations. 
As discussed in the previous section, for a mean velocity of 
~3 km s _1 , the depletion will be ~10 5 years. 

The energy of any kind of turbulence will decay if there is 
not any source injecting new energy. It is interesting to compare 
the depletion time scales with the turbulence decay time scales 
to establish if the turbulence can maintain mean velocities 
that are high enough to match the depletion ti me scales. The 
turbu l ence decays as r g , where a~0. 8-1.0 jMac Low et alJ 
Il998l IStone et alJll998l) . The decay time scale of these tur- 
bulences is similar to one dynamical time (td yn =L/cr v , where 
L is the turbulent length and q- v the turbulent velocity, 
Avila-Reese & Vazauez-SemadenilbOOlt lElmeereen & S calo 
2004 . br both HD and MHD turbulence. The dynamical time 
for turbulence in the GC, assuming a turbulent length of 5 pc 
and a turbulent velocity dispersion of 15 km s _1 (observed GC 
linewidths), is ~3xl0 5 years, close to the depletion time scale. 
Since we do not observe large depletions in our data, the ex- 
istence of faster shocks with velocities high enough to pro- 
duce the sputtering of the grain mantles and with the time 
scales between shocks similar to the time scales for depletion, 
~10 5 years is required. Shocks with velocities of ^;6 km s are 
required in order to sputter molecules from the icy mantles on 
graph ite and silicate grains by heavy atoms JPraine & SalpeteJ 
I1979I) . Then, if due to a sudden event, the turbulence generated 



by this event could explain the presence of complex organic 
molecules in the GC clouds for ~10 5 years. 

It is so far unclear what produces the supersonic turbulence 
and/or frequent shocks in the GC. The high-velocity shocks 
and the presence of a highly turbulent medium in the GC could 
be due to: a) the kinematics of the gas subjected to the barred 
potential of our Galaxy, b) the wind-blown bubbles produced 
in evolved massive stars like supernovas, or c) cloud-cloud 
collisions. The time scales needed in the GC to maintain the 
large abundances are much smaller than those associated to 
the quasi-circular orbits (~lxl0 7 years), and a large number 
of large scale shocks would be required to explain the ob- 
served properties of the complex molecules. The other possi- 
ble mechanisms could explain the observations if cloud-cloud 
collision and energetic events driven by massive stars occur 
on time scales of 10 s years. In the scenario of a recent star 
burst in the GC as proposed to explai n the fine structure line s 
of ionized gas (Rodrig uez-Fernandez & Martin-Pintadol20 05i. 
frequent energetic events associated to massive stars are ex- 
pect ed to produce shocks with moderate velocities of ~10km 
s" 1 iMartin-Pintad oet al J 19991) ejecting complex molecules to 
gas phase and a large amount of energy into the ISM. 

6.4. The chemistry of complex molecules in the 
Galaxy 

The proposed gas-phase chemistry for some complex 
molecules does not seem to explain the similar abundances ra- 
tios observed in the GC clouds, in hot cores, and in hot corinos. 
From the empirical data gathered there, we can conclude that 
the most likely explanation for the large abundance of CH3OH, 
HCOOH, HCOOCH3, (CH 3 ) 2 and C 2 H 5 OH is that all these 
molecules are ejected from grains. It is interesting to note that 
in this scenario, the ice-mantle composition in complex organic 
molecules in the GC and in hot cores must be very similar. 
This is surprising in view of the changes of the complex or- 
ganic molecules abundances observed in the hot corinos. The 
presence of complex molecules on grains can be due to grain 
chemistry or to the depletion after formation in gas phase. 

Further investigations of the abundance of organic 
molecules in dark clouds are needed to constrain the molecules 
formed on ice mantles by grain chemistry in different types of 
molecular clouds. 

6.5. Sources with low C 2 H 5 OH abundances 

We also observed some clouds where the relative abun- 
dance of some organic complex molecules like C2H5OH 
seems to be significantly smaller than expected from 
these observed in other sources in the GC. These sources 
are the Sickle (MC G+0.18+0.04 and MC G+0.20+0.03) 
and the -30 km s _1 velocity component of the TRA 
(MCG+0. 13+0.02). We detected high CH 3 OH abundances 
of ~10~ 7 , but only significant upper limits, ^0.4-5 xl0~ 9 to 
the C2HgOH abundance. The stars from the Quintuplet and 
the Arches Clu sters are said to heat and i onize t he Sickle 
and the TRA ( Rodrig uez-Fernandez et alJ bOOlbl) creating 
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photodissociation-regions (PDRs) around them. In particular, 
the -30 km s" 1 com ponent is clearly associated with the TRA 
JSerabvn & Gustenll987l) . 

The chemistry in these sources is then expected to be 
strongly affected by the presence of the UV radiation. It is very 
likely that the complex molecules in these sources are ejected to 
gas phase by shocks with the same abundances as observed in 
other GC molecular cloud. Additional evaporation from grain 
mantles can be ruled out since the heating will require UV pho- 
tons that will also destroy the complex molecules. Due to a 
large UV radiation-field, the relative abundance of the com- 
plex molecules will be affected by dissociation. One can make a 
rough estimate of the effects of the UV radiation on the molec- 
ular abundances by comparing their photodissociation rates. 
Using the UMIST RATE99 dLeTeuff et all2000l) and the Ohio 
State University databases 1 , we found that the photodissocia- 
tion rates of (CH^^O and C2H5OH are larger than those of 
CH3OH, HCOOH and H 2 CO. 

In Fig. |9]we show the results of a simple model in which 
we represent how the ejected complex organic molecules are 
photodestroyed as a function of time and visual extinction. The 
model considers a far-ultraviolet incident field 10 3 times higher 
than in the local ISM (Gq=1 3 ), which is the typical radi ation 
field observed in GC PDRs (Rodriguez-Fernandez et al. 2004). 
For the regions with low A v , all the abundances decrease, but 
the abundance of CaHgOH decreases faster than the other or- 
ganic complex molecules. Obviously the time scale for pho- 
todissociation depends on the visual extinction. We estimate a 
visual extinction of A v ~4-17 mag for the Sickle and the TRA 
using our H2 c olumn densities in Tableland the conversion 
factor given bv iBohlinl dl975l) . A v (mag)=10" 21 xN H2 (cm- 2 ). 
The photodissociation times scales for CH3OH will be compa- 
rable to depletion time scales for an averaged visual extinction 
of ~8 mag. 

For this visual extinction the other complex organic 
molecules will photodissociate much faster than the depletion 
time scale. This suggests that in the presence of UV radia- 
tion, large abundances of molecules like C2H5OH would re- 
quire more frequent shocks than in the shielded regions. This 
is illustrated in Fig. [H)] where we present the results of a 
simple model to show the time evolution of complex organic 
molecules in a region of A v =8 mag. In this model, we consider 
the effects of the UV radiation on the complex molecules that 
have been ejected by periodic shocks of 5xl0 4 years. They are 
ejected with the abundances measured in the GC and are de- 
stroyed because of the photodissociation. This simple model 
reproduces the constant high abundances observed for CH3OH, 
while the abundances of C2H5OH can vary by ~3 orders of 
magnitude. In agreement with the observations, the photodisso- 
ciation can reproduce the low abundances of C2H5OH in spe- 
cific regions in the GC clouds. Further observation of complex 
organic molecules toward PDRs in the Galactic disk and in the 
Galactic center are required to test the proposed scenario. 



1 www.physics.ohio-state.edu/~eric/research.html 
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Fig. 9. Photodestruction of different complex molecules as a 
function of time for different visual extinctions (A v ). The initial 
molecular abundances for the molecules are considered to be 
those typicall y measured in the observed GC c louds. Go=10 3 , 
observed by Rodri guez-Fernandez et alJ d2004l) . is used for the 
UV radiation field. 

7. Conclusion 

We have found very high abundances for the complex or- 
ganic molecules CH 3 OH, C 2 H 5 OH, HCOOCH3, (CH 3 ) 2 0, and 
HCOOH in the GC regions, similar to or even higher than those 
measured in Galactic hot cores. The abundance ratios of these 
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Fig. 10. Time evolution of the abundance (top) and relative 
abundance (bottom) of complex organic molecules with respect 
to that of CH3OH in a PDR with A v =8 mag where a periodic 
shock every 5 x 10 4 years enhances the abundances to their ini- 
tial values. The radiation field is Gn=10 3 . 

molecules are fairly constant in all the GC clouds. The most 
likely explanation for these large abundances is that they have 
been ejected from grain mantles by shocks. The highly turbu- 
lent molecular clouds in the GC could help to maintain the large 
gas phase abundance of complex molecules for time scales of 
~10 5 years after the ejection by shocks. The rather uniform 
abundance ratios in clouds distributed over 200 pc indicate that 
the average composition of grain mantles is similar for all the 
GC clouds and also for the hot cores in the Galactic disk. We 
found that, in the GC PDRs, UV photons can differentially 
photodissociate the complex molecules, changing their relative 
abundances dramatically. 
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Table 6. Relative abundances and H9 column densities 
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1 Derived from CH 3 OH. 

2 Derived fro m C 2 H 5 OH. 

3 From 13 CS jMartm et all2006l) . 

4 Derived from HCOOCH3. 

5 Assumed 8 K. 

6 Derived from (CH 3 ) 2 0. 
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